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ABSTRACT 

We have made a detailed 21 cm study of areas that have the smallest known 
amount of HI In the northern sky. These observations were corrected for stray 
radiation using a method described in an Appendix, and have an estimated 
uncertainty in Ngi of <5 x 10^® cm”2. The region of main Interest, around 
a ■ 10^45®, 6 ■ 57“20' , has a minimum Njjj of 4.5 x 10^® cm“2. Spectra taken at 
21* resolution over a field 4* x 3" In this direction show up to four HI line 
components. TVo, near 0 and -50 km s“^, are ubiquitous; their combined column 
density Is always >4.5 x 10^^ cm"2. There Is also a narrow component at 
-10 km s"^ attributable to a diffuse cloud covering half of the field, and 
scattered patches of HI at v < -100 km s“^. The low and intermediate velocity 
components have a broad line width and are so smoothly distributed across the 
region that it is unlikely that they contain significant unresolved angular 
structure. Eight other low column density directions were also observed. 

Their spectra typically have several components, but the total column density 
Is always > 7 x 10^^ cm“2 and changes smoothly along a 2* strip. Half of the 
directions show narrow lines arising from weak diffuse HI clouds that contain 
0. 5-3.0 X lOl® cm-2. 

Directions north of declination -40* that have Ngj < 7 x 10^® cm“2 appear 
to be quite rare, and it is unlikely that the HI column density ever drops 
significantly below 4.5 x 10^® cm"2 anywhere in the sky. Thus the ISM provides 
an effective barrier to the direct passage of ionizing radiation at E < 100 eV, 
either from extragalactlc sources Into the galactic plane, or from the solar 
neighborhood Into the galactic halo. 


3 


I. INTRODUCTION 

This paper contains the first results of a study of HI emission In the 
regions of the northern sky that have the least neutral hydrogen. These areas 
are of particular Interest because they often have so little HI that spectra 
can be decomposed, v^th confidence. Into a few Individual features that 
persist over degrees. This Is very different from the situation In directions 
where HI Is plentiful and spectra are obviously a confused blend. We hope 
that the low column density regions will provide a starting point for 
understanding the small-scale organization of the Interstellar medium. (Xir 
second goal, and one that Is stressed In this paper. Is to search for 
directions that might have extremely small amounts of neutral hydrogen and 
thus provide transparent lines of sight to extragalactlc objects at EUV and 
soft X-ray wavelengths. 

Low HI column density directions were chosen guided by the HI survey of 
Stark et al . (1986; hereafter SBLH), which was made with the Crawford Hill 
hom-reflector of the AT&T Bell Laboratories. This survey has an angular 
resolution of only 3* x 2” , but It Is uncontaminated by stray radiation and 
thus able to supply large-scale HI column densities more accurately than 
previous studies. Our primary target was the region around t > 150*, b > 53*. 
It has the lowest column density In the SBLH survey, and, because It appears 
as a local hole In column density maps, we refer to that area as "The Hole." 
Other directions studied Include 1 ■ 86*, b ■ 76*, the second lowest Nm 
region In the SBLH survey, the north and south galactic poles, and several 
regions which were found to have a low total N^x In the course of observing 
> 100 high latitude directions (Jahoda et al . 1985; Lockman, Hobbs and Shull 
1985). 
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The areas that were observed are listed in Table 1 with their over 

the SBLH beam. Our higher resolution 21' observations have revealed a few 
directions with a total %x of only 4.5 x 10^^ cm"2, the lowest ever reported 
in galactic 21 cm studies, but in general, directions with N^x < 7 x lO^^ cm"2 
seem quite rare. 

Accurate measurement of low brightness galactic HI profiles requires that 
stray radiation be removed from the spectra. The technique used for this, and 
a discussion of ways to Improve the correction for instrumental baselines, are 
given in Appendices. 


1 


5 

II. OBSERVATIONS AND DATA REDUCTION 

All of our 21 cm data were obtained with the 43 m (140 foot) telescope of 
the NRAO In Green Bank. West Virginia, during April 1984. At the frequency of 
the 21 cm line this telescope has a half-power beam-width (H?BW) of 21'. The 
receiving system consisted of dual cooled FET amplifiers and a hybrid mode 
feed that gave simultaneous linear polarizations (Norrod 1984). The total 
system temperature, on the telescope at zenith, was 23 K. the aperture 
efficiency was 60Z and the main beam efficiency was about 801. Spectra 
covered 1 240 km s~^ about zero velocity LSR at a resolution of 1.0 km s~^. 

A parabolic baseline was removed from each spectrum using a technique 
described In Appendix B. The noise In the final spectra Is typically 
< 0.07 K (T^) and a typical error In a column density determination. Including 
baseline uncertainties, is 5 x 10^® cm”2. Column densities were derived under 
the assumption that the lines are optically thin. This Is a good assumption, 
as the peak line temperatures are always < 3.5 K (Tj)) and typically < 2.5 K. 

One-dlmenslonal maps of the HI distribution were made by measuring spectra 
every 10' In right ascension along 2* strips at constant declination. Toward 
The Hole a number of strips were combined to give complete coverage of an area 
4* X 3* at 10' X 10' spacing. 

The treatment of stray radiation Is discussed In detail In Appendix A. In 
brief, an area Including the entire main beam of the Crawford Hill 
horn-reflector of AT&T Bell Labs was mapped with the 43 m telescope. This map 
was convolved with the horn-reflector beam pattern, and differences between 
the convolved 43 re spectra and that actually observed using the Bell Labs 
antenna (SBLH) are Interpreted as arising from stray radiation which la then 
removed from Individual 43 m spectra. The stray radiation correction was 
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redetermined at least every 1*5 hours and independently for every direction. 
We are confident that the accuracy of the resultant spectra is limited by 
noise and uncertainty in the instnimental baseline rather than by stray 
radiation. 
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III. THE HOLE 

Figure 1 shows the low column density end of the SBLH HI survey binned by 
the total column density In each spectrum. Of the 22000 spectra In that 
survey, which covers 6 > -40*, only 0.9% have a total HI column density 
< 1 X 1020 cm“2, and directions with Nm < 7 x lO^^ cm”2 comprise < O.IZ of 
the survey. These lowest column densities (shaded portion of the figure) all 
come from a part of the sky centered roughly on o - 10|?7, 6 ■ 58*. We call 
this region "The Hole." 
a) The Total Column Penslty 

A map of the Integrated HI column density over •'12 square degrees of The 
Hole is shown In Figure 2. The map is centered at 01959 ■ 10^ 45®, 

61950 * 20'. Within the area of complete coverage (Indicated by the 

dashed border) spectra were taken every 10' In both coordinates. Each 
position was observed at least twice; a total of about 1300 spectra were 
combined to make the map. The contour levels are In units of 10^^ cm"2 and 
the minimum Njn, 4.5 x 10^^ cm”2^ Is observed at offset position -60' , O' and 
again near +75', +50'. In retrospect, the map should have been extended to 
lower right ascension, but the observed spectra were all heavily contaminated 
by stray radiation (see Appendix A) so there was no way to know what the true 
signal was during the observations. Nonetheless, we have a few additional 
spectra outside the borders of Figure 2 which Indicate that the main minimum 
In Ngi Is covered In our survey. Therefore, in all likelihood. Figure 2 shows 
the smallest HI column density directions north of declination -40*. 

The HI spectra toward The Hole are remarkably simple. Examples are given 
In Figure 3 Identified by their offset positions for reference to the map In 
Figure 2. Each spectrum can be decomposed Into at roost four well-defined 
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Gaussian components, and Figure 2 was derived from summing their Individual 
contributions. The component with the largest column density Is broad and 
near zero velocity. Every spectrum also contains another broad line at an 
Intermediate velocity near -50 km s**^. These two components contain most of 
the column density In The Hole. In some directions the spectra show a third 
component that is a bright narrow line attributable to a diffuse HI cloud. It 
covers about half of the region. The fourth, high velocity, coraponant shows 
up as occasional weak emission at v < -100 km s~^ that Is never brighter than 
the example shown In Figure 3. 

The spectrum from offset (-65*, O') shown In Figure 3 Is the weakest 
galactic HI line detected to date, both In terms of peak brightness 
temperature (0.6 K T^,) and total Ngi (4.5 x 10^^ cm“2). Nonetheless It 
contains two distinct components and both have a line width (FWHM) '^25 km s~^. 
b) The Individual Components 

Column density maps for the spectral components are given In Figure 4. 

The low velocity emission is the most uniform spectral feature over The Hole. 
It contains half or more of the total Ngj nearly everywhere Inside of the 
mapped region. Its column density never falls significantly below 
3 X 1019 cm~2 even over the broad minimum that runs diagonally across the map, 
and it reaches a maximum of just over 4.5 x 10l9 cm“2 at the edges of the map. 
Its peak brightness temperature ranges between 0.5 and 0.8 K, while Its width 
(FWHM) is always close to 30 km s~l and has no significant trend across the 
map. The column density structure thus reflects mainly changes In the peak 
temperature. Over most of the map the peak velocity Is close to -2 km s“l. 
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There is c trend to lower velocities, but by only s few ka s**^, in a rou;;hly 
triangular area whose apex la near offset 0,0 and whose base extends over most 
of the bottom of the map* 

The Intermediate velocity component ia alao ubiquitous. Its minimum 
column density is only slightly less than 1 x 10^^ cm'^ while its maximum is 
slightly greater than 3 x 10^^ cm"2 at the map edge. The average line width 
is about 30 km s"^, but it decreases roughly with right ascension from about 
35 km s~^ at the right edge of the map to about 22 km at the left edge. 
Most of the column density structure is caused by changes in the peak line 
temperature which ranges between 0.2 and 0.7 K. In a few directiona this 
intermediate velocity feature is as bright e.s the low velocity feature, but it 
always has less column density* Ihe peak velocity varies from about 
-45 km s~^ at low right ascension to about -51 km s~^ at high right ascension 
but with considerable scatter within these limits. The velocity gradient is 
opposite to the sense of galactic rotation. 

The location of the cloud is shown 1:* Figure 4, and maps of its line width 
and peak velocity are shown in Figure 5* Column density contours are in 
10^^ cm~2, and the line width and velocity contours are in km s"^* The dashed 
contour shows the extent of the minimum detectable emission* The cloud 
consists of northern and southern parts that are Joined by faint emission* 

The northern part (A6 > -30* ) has a narrow line width: it Is usually <6 km 

s**^ and occasionally drops slightly below 4 km s~^. There is also a gradient 
in line width across the northern part: Av generally increases with right 

ascension and is larger at the cloud edge than near the center* The line 
width of the southern component is much larger with a stronger gradient* The 
width is always > 6 km s~^ and rises to Just over 14 km s~^* The two cloud 
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conponentf alto differ in their velocity. The northern component hat a alight 
eaat-weat gradient of a few ka a**^. while the velocity of the southern 
component decreasea rapidly to the southeast. It is possible that the 
southern component is actually a superposition of two sirall clouds. 

Even though the column density in the cloud peaks in the areas where the 
main low velocity feature is also strongest, there is no clear connection 
between these two spectral components. The velocity end line width gradients 
in the cloud have no counterpart in the low velocity emission, and the cloud 
always has a more negative velocity by up to 10 km a~^. 

The cloud is much fainter chan typical diffuse HI clouds which have column 
densities x 10^^ cm~^ (e.g.. Spitzer 1978). but is very much like the 
"cloudlets'* discovered by Heiles (1967). If it is '^100 pc away, then it is 
only a few pc across and has a mean density *^1 cm~^. 

High velocity emission over The Hole is weak and scattered. Its column 
density (Fig. 4) is never much more than 1 x 10^^ cm~^. Dashed contours mark 
the maximum but uncertain extent of the high velocity patches. Most of the 
gas has a velocity of about -105 km s~^ except in the northeast corner of the 
map where the velocity is near -130 km s~^. Line widths are about 25 km a~^ 
and the peak brightness temperature ia never greater than about 0.2 K. Qouda 
at similar high velocities have been reported in this general part of the sky 
(e.g.. Giovanelll 1980). 
c) Spatial Structure in the Emission 

The dominant HI components toward The Hole look very smooth in the column 
density maps. A quantitative measure of smoothness is the first order 
structure function D: 

D(6) 5 i J [Nhi(1) - Hhi (l+0)]2 
" n 


( 1 ) 
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where 6 le en angular separation and the sun Is over the n sppropriate spectra 
In the map (see, e.g., Rutman 1978; Simonettl, Cordes, and Heeschen 1985). 

The first order structure function gives the mean square difference In 
over the entire map, of points separated by the angle 6. Because our map la 
sampled every 10* while the 43 m beam la 21* , the structure function for 
6 ■ 10' should contain almost entirely noise and baseline uncertainties which 
are unrelated to the spatial emission structure, l.e*, D(IO') * 2o^. The 
uncertainty In a total column density, estimated this way, la 
o - 0.4 X 10l9 cm~2, which Is about what Is expected from our error analysis 
(SII). It thus makes sense to define a reduced structure function 
D'(6) = [d( 8) -D(10')]^/^ which contains mainly Information on the spatial 
variations In Nhx* Figure 6 shows D* for the total N»t and for the dominant 
low and Intermediate velocity components, noxTiallted by the typical column 
density of each: 5.5 x 10^^ cm"2, 3.5 x lO^^ cm”^, and 2 x 10*^ cra“2, 

respectively. This figure thus gives the proportional average rms difference 
In column density over a given angular scale. It shows that nearby points are 
highly correlated, distant points less so. 

The fact that D* (6) Is nearly monoLonlc confirms what the appearance of 
the column denalty maps suggests — that the column density changes are 
dominated by gradients, not by relatively small scale Irregularities. The 
apparent detail In the total column density map Is somewhat misleading In thla 
regard, for what Is not readily apparent there, but shows up In the structure 
function analysis. Is that there la always an underlying constant Ngx; 
the maximum difference between any two points In the map Is always less than a 
factor of 2, and la usually less than a factor of 1.25. Although the 
Intermediate velocity component can show fluctuations of up to 40Z at large 
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■calet* th« low velocity, and dominant, component la quite smooth on small 
scales. For these reasons, and because most of the column density is in 
components that have a line width <^25 ka s~^, %re think that significant 
structure in the total Ngj angular scales smaller than those observed here 
is unlikely. This issue is discussed further in the following section. 
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IV. OTHER LOW Nhi DIRECTIONS 

Neutral hydrogen obaervatlona, corrected fur stray radiation, were made In 
right ascension strips 2* wide across eight areas of low Nj^j In addition to 
The Hole. Figure 7 shows the distribution of Nhj for these strips, snd 
Figure 8 ^Ives the spectrum at the center of each strip. One contaon 
characteristic of every field Is that the emission varies sm.'othly and has no 
local peaks or holes. In the example of greatest emission variation, the 
strip centered at 1 - 182*, b ■ 4-58?4, much of the change Is caused by 
structure In a well known high velocity cloud that lies In this direction 
(e.g., Glovanelll 1930). The dashed line In the Ngj curve for this strip 
shows the hydrogen omitting the contribution from the HVC. None of the eight 
areas ever has an Ngx < 7 x 10^^ cm''^. 

The spectra always contain one or more broad components and often a 
relatively narrow, bright line like that from the cloi* * which covers part of 
The Hole. These “clouds" have column densities In the range 0.5 to 
3.0 X 10^^ cm”2, llnewldths (FWHM) of 4-12 km s"l and velocities between -15 
and -fl km s~^. They cover about half of the low f'gx directions. These 
relatively «reak features, which would be difficult to disentangle from the 
bright HI emlsslcs In more typical hlgh-latltuda directions, seem to pervade 
the Interstellar medium. As In The Hole, these clouds are not obvloualy 
connected with the main low velocity gas. 

The general HI atnicture In regions of low Ngx can be described by several 
statistics. For each strip we have calculated the mean Ngx >nd Its dispersion 
0. The ratio o/ <Ngx> ranges between 2.3Z and 11. 5X with an average value of 
7.2Z. The reduced and normalized structure function, 0'(6), as described In 
fllle, averages 0.10 ± 0.06 at 6 - 60* and 0.15 ± 0.09 at 0 - 90*. Valuas of 
D* for the total Ngx In the Hole (Fig. 6) are within the lo range of these 
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results. Some directions are remarkably smooth In total The strips 

centered at 5620 + 66tl and 85*6 + 7529 have o/<Njjj> of 2.3Z and 3.7%, 
respectively; toward 5620 + 6621, where the total Nhj Is always 
< 1 X 10^0 cm“2, the reduced structure function D'<60’> - D'<90'> - 0.02. 

Even the cases of most extreme spatial variation, the south galactic pole and 
18220 + 5824, have a a/<Nm> that Is only 11. 5Z and 11. It, respectively. 

Another perspective on the essential smoothness of in these areas 
comes from comparison of over the ^0.7 square degrees of our strips 

with the. average over the 6 square degrees of the SBLH survey beam 
(Table 1). For the eight areas the average ratio of our column densities to 
those of SBLH Is 1.01 ± 0.09 (la). The maxltn»im difference In column density 
between a single spectrum at 21’ resolution and one at ~3® x 2* resolution Is 
5.8 X 10^^ cm“2, and this Is caused by a small high velocity cloud In the 21' 
spectrum. If high velocity gas Is neglected, then the maximum difference Is 
half this value, only 3 x 10^^ cm”^. In low Njji directions only HVC's are 
significantly clumped on angular scales <2*, and coarse resolution surveys 
provide a good guide to the total column densities seen at higher resolution, 
within a la uncertainty of '‘lOZ. A similar strong correlation In between 
high latitude positions separated by 8' has been observed by Payne, Salpeter, 
and Terzlan (1983). 

These considerations suggest that the hydrogen Is extended and not 
contained In very small, unresolved, clouds. For example, there would have to 
be 100 small clouds In the 0.2 square degrees of the 43 m beam to produce 
the observed a/<N{{x^ 0.1. At a distance of 100 pc, each would have a 

diameter 0.1 pc and a density <n> *>■ 200 cm"^. 
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V. SUMMARY 

Areas of low total neutral hydrogen at 6 > -40* have a number of common 
features: 

1. Directions with < 7 x 10^^ ca"2 are rare. We have found only one 
such area: The Hole near t 150*, b ■ 53*. Even there the column density is 
> 5 X 10^^ cm“2 except over two small regions where it reaches 

4.5 X 10^^ cm"2, 

2. The total column density, even in these low regions, varies quite 

smoothly with position. A typical rms dispersion in ^or all points in a 

2* long strip, is only 101 of the mean. The rms difference between column 
densities 1* apart on the sky is < 15Z of the mean. In regions of low Nyj the 
smallest angular scale features are high velocity clouds. These contribute 
most of the small-scale structure in some directions. The part of The Hole at 
offset (-45*, -50') in Figure 2 is one such example, as is the area around 
182!0 + 58J4 (Fig. 8). Large changes in Ngx with position typically arise 
from broad gradients, as emphasized by Payne, Salpeter and Terzlan (1983), and 
not from small-scale fluctuations. Thus observations at the rather coarse 
angular resolution of ">3* x 2* (Stark et al. 1986) give a good estimate of the 
total HI which would be seen at higher resolution. 

3. Profile shapes in low %x regions are usually complex. Most show 
between three and five components over a few degrees of the sky. But most of 
the emission is in broad spectral features with a FWHM > 25 km s~^ that change 
smoothly with position. 

4. About half of the low N^x directions, including The Hole, show a 
relatively narrow spectral component with a FWHM between 4 and 12 km s'^, like 
that expected from a small diffuse cloud. The features contain 

0.5-3.0 X 10l9 cm~^ and have velocities -15 < v < 1 km s~^. They do not seem 
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to be simply connected with the main, broad, lev velocity emission lines. 

These %roak, narrow features would be difficult to detect In the brighter and 
more confused HI profiles that are typical of most of the sky. They are 
probably the faint end of the ’'cloudlet” population discovered by Helles 
(1967). 

Although our observations are limited to 5 > -40* , It Is likely that they 
cover the lowest Ngx directions anywhere. The all-sky HI data of Cleary, 
Helles and Haslam (1979), which contain varying amounts of stray radiation, 
still show a number of directions at 6 > -40* with Ngj < 1 x 10^® cm”^^ but no 
such directions south cf -40*. It therefore seems probable that column 
densities significantly less than our minimum, 4.5 x 10^^ cm~^, do not exist 
In any direction, and that the Interstellar medium Is opaque (t > 2) to 
extragalactlc radiation at energies between 14 and 100 eV. 

This research was supported In part by the National Aeronautics and Space 
Administration under grant NAG 5-629. We are grateful for the use of the 
Midwest Astronomical Data Reduction and Analysis Facility (MADRAF). 


17 


APPENDIX A 

CORRECTION OF 21 cn SPECTRA FOR STRAY RADIATION 

Radio telescopes have weak sldelobes that cover the entire sky. Neutral 
hydrogen emission entering these sldelobes adds to emission detected via the 
main beam. When the main beam lies on an area of weak HI while the 
sldelobes cover bright emission from the galactic plane, the ‘’stray" radiation 
can comprise half of the observed signal. This problem has been recognized for 
some time (e.g., Van Hoerden, Takakubo and Braes 1962), and elaborate methods 
have been devised to reduce or remove Its Influence (e.g., Kalberla, Mebold, 
and Reich 1980; hereafter KMR). 

Figure 9 shows an example of the problem. Ihese two spectra were obtained 
with the NRAO 43 m telescope at Green Bank from Identical observations of 
Identical parts of the sky at two different hour angles. The 43 m telescope Is 
equltorlally mounted, so when It tracks an object Its sldelobe pattern Is fixed 
with respect to the sky. Even so, the stray, radiation changes as far sldelobes 
move above and below the horizon. This, together with the fact that stray 
radiation enters a spectrum at a velocity which Is not reduced to the LSR, 
means that the stray radiation spectrum Is a function of telescope hour angle 
and elevation, and the time of year. But, as emphasized by KMR, the time 
varying component of stray radiation, which Is so prominent In Figure 9, Is 
only a lower limit to the total amount. 

Stray radiation can be removed from a spectrum through an all-sky 
deconvolution provldad that the complete antenna response pattern Is known to 
below the -50 dB level. This method was used by KMR, who modelled the far 
sldelobes of the 100 m Effelsberg telescope to the -70 dB level (the far 
sldelobes of most antennas cannot. In general, be measured accurately) and 
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estimated a stray radiation spectrum by convolving the antenna pattern with 
the *> 2 X 10^ spectra in the Weaver and Williams (1973) and Heiles and Habing 
(197A) HI surveys. We found this technique unattractive primarily because the 
far sidelobe pattern is sensitive to the illumination pattern and geometry of 
a telescope. Any change in the receiver feeds, or the telescope structure, 
would require a laborious recalculation of the sidelobe response. Such 
changes, particularly in the form of receiver improvements, occur routinely at 
Green Bank. 

Some radio telescopes, such as horns and hom-ref lectors, have negligible 
far sidelobes and hence little or no stray radiation. One such instrument is 
the Crawford Hill horn-reflector of the AT&T Bell Laboratories (hereafter the 
BL telescope) which was used recently for an HI survey of the northern sky 
(Stark et al . 1986; hereafter SBLH). Although this survey has minimal stray 
radiation, the telescope has an angular resolution of only >^3* x 2° at 21 cm; 
for many purposes HI needs to be measured on smaller angular scales. 

We have developed a method of "bootstrapping'* A3 m HI spectra to the SBLH 
survey which yields data at 21' resolution uncontaminated by stray radiation. 
The method has been used for several years with receivers having quite 
different illumination patterns, and it gives repeatable, reliable results. 

This Appendix contains all the information that is necessary to use the 
method on any telescope. Section a) describes the technique and the 
fundamental assumptions on which it rests. Section b) contains a discussion 
of the BL antenna and of the assumption that the SBLH survey data have no 
significant stray radiation. In section c) we discuss the methodology of the 
observations, and show, in d), the typical levels of stray radiation that are 
present in A3 m spectra. Extensive tests of the other basic assumptions are 
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reported in section e) and in f) we summarize the technique and discuss its 

limitations. 

a) The Technique 

l.t principle, a high angular-resolution map made with a large 
fllled-aperture telescope (the NRAO 43 m) can be convolved with the beam 
pattern of a smaller telescope (the BL horn-reflector) to reproduce the 
low-resolution observations exactly. We use this fact to interpret any 
differences between a convolved 43 m spectrum and the one observed at BL as 
originating in stray radiation in the 43 m data. Helles (1974) made similar 
calculations in a few directions to determine a rough stray radiation 
correction to the Heiles and Habing (1974) HI survey. 

The stray radiation spectrum, <S(v, oq, averaged over '^3“ x 2* 

toward a particular direction at a particular time, is defined as 

<S(v,oo,5o)> - <TnraoCv,cxo»5o)> " ^bl (v,ao,6o)» A1 

where 

<TnrAo(v»“O»®0)> ■ / TnrA0Cv»®»*) f(a-oto,fi-5o) A2 

Q 

the are spectra from a map made at the 43 m, TBj^(v,aQ,6Q) is 

the spectrum observed by SBLR, and f(a-aQ,6-6Q) is the antenna pattern of the 
BL hom-reflector normalized so that 

/ f(a-oo,fi-«o) dn - 1. A3 

Q 

The stray radiation spectrum is not only a spatial average over the BL beam, 
but an average over the time. At, required to obtain the 43-m spectra. If we 
assume that the stray radiation contribution is essentially constant over some 
area around (o(),6q), and that it is constant over At, then the average stray 
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radiation spectrum derived from a A3 m map can be removed from Individual 
spectra In that map to produce "corrected," hlgh-resolutlon spectra: 
Tcor(v.o,5) - TnRA0(v,a,6) - <S(v,ao,6o)>* 

The accuracy of depends on the validity of three assumptions: 

(1) That the BL telescope spectra are correct and contain no significant 

stray radiation. Note that all differences between and are 

transferred to T^^^. Those Include noise, the Instrumental baseline, and 
the relative calibration of the spectra. 

(2) That the A3 m stray radiation spectrum has not changed significantly over 
the time. At, required to make the A3 m map. 

(3) That the A3 m stray radiation spectrum does not change significantly over 
the area of the BL main beam, i.e. that S(v,a,6) ■ <S(v,ao,5o)^ 

Io-OqI < 1“ and |A-5ol ^ 

Tests of the assumptions are described In detail in sections b) and e). 
b) The BL Antenna and the SBLH Survey 

The BL hom-reflector at Crawford Hill is described by Crawford, Hogg, and 
Hunt (1961). Its beam shape can be calculated from its far-fleld pattern, 
which has a directional dependence given by 


e(0l»®t) 


i/ dC / d^ Eg(^,C) «•> exp[-l(Dic (cos Z sIoBl + sin C sin 
”^o “♦o 


0t)] 

AS 


where u - cos ^/(1-sln^), k - hutfXt ^ ^s the wavelength and f, the focal 
length of the reflector. Is 20 feet (6.096 m). The Integration Is over the 
longitudinal and transverse r.ngles within the horn, ^ and respectively, 
between the half flare angles * 1AS036. Angles In the far field are 
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from the beam center In the plane containing the horn axis ( 6 |^) and the plane 
perpendicular to the horn axis ( 67 ). For linear polarization the vector 

four components, two In the longitudinal direction 


Q z) 


cos C sin y cos 
and two In the transverse direction 

sin 5 sin y + cos 




Ot) 

OU 


sin C cos y A 6 


cos C cos y 


A7 


where y Is the angle between the polarization In the horn and the ^ • 0 plane 
and Is zero for purely transverse polarization. 

The beam shape of the hom-reflector Is a function of the antenna 
polarization and the HFBW varies from about 2?9 x 1?9 In transverse 
polarization to 2?2 x 2JA In longitudinal polarization. The polarization 
angle of the SBLH survey varied systematically with elevation, being purely 
transverse at an elevation of 45”. Thus the antenna pattern must be 
calculated separately for each declination. 

Almost all SBLH spectra were taken at transit; observations north of the 
zenith were usually made at an azimuth of 0 ” , south of the zenith at an 
azimuth of 180*. Thus 6 ^ and By correspond to offsets In apparent a and 6 , 
respectively, and the BL antenna pattern, f(a-ao, 8-5q), can be written 


f(o-oo» ■ const, |E(a-oo, S-Bq) P 


A8 


where the constant provides the normalization (eq. A2). 

The calculated pattern does not deviate from the observed horn-reflector 
pattern (Crawford, Hogg, and Hunt 1961; Penzlas, Wilson, and Encrenaz 1970) 
by more than a few percent, except that In longitudinal polarization 
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two far sidelobes, arising froir spillover, and from diffraction over the drive 
wheel, appear at a level of up to -35 dB. The area under these far sidelobes 
Is < O.iZ of that In the main beam. We estimate that in the worst case these 
far sidelobes could contribute as much as 8 x 10^^ cm~^ to the observed 
column density, but a more typical value would be 0.5 x 10^^ cm'^. 

More than 99% of the BL antenna response lies within 14* of the main beam, 
which Itself contains 92Z of the tot-'' response (Penzlas, Wilson, and Encrenaz 
1970). It would be a reasonably straightforward task to "correct” the SBLH 
survey for Its own near sidelobes, but. It appears that these corrections 
would be Insignificant In nearly all directions. For example, toward the 
minimum column density position In The Hole the correction for sidelobes 
within 14* of the main beam would be less than 2 x 10^® cm”^, or about 4Z of 
the uncorrected column density. This Is an extreme case, because here the 
column density In the main beam Is very low and all sidelobes lie on higher 
Nhi regions. The near sldelobe contribution to the SBLH survey' thus seems 
insignificant. 

The antenna response more than 14* away from the main beam, except for the 
unique spillover and diffraction lobes. Is difficult to locate experimentally. 
Our calculations Indicate that there should be two long, weak, sidelobes, 
extending east-west In longitudinal polarization and north-south In transverse 
polarization. If these lobes contained as much as 0.5% of the total response, 
we estimate that they could typically contribute 2 x 10 cm"2 to the observed 
column density at high latitudes. 

It appears that the stray radiation component of the SBLH survey should 
generally be less than a few times 10^^ cm'^, except for Isolated Instances 
where It might reach 1 x 10^^ cm"^. Thus, to a relatively high degree 
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of accuracy, we can treat the SBLH aurvey as If It had no sldelobes, and 
perform the convolution of the 43 m data (eq. A2) only over the area of the 
main BL beam. 

c) Making the Observatiom 

A 43 m map made over an elliptical area 5t3 x 4** always covers the main BL 
beam to below the lOZ response level and provides an adequate estimate of the 
BL spectrum. Moreover, a complete sample of this area at the 21' resolution 
of the 43 ra, which would require 550 spectra, is not even necessary. We have 
found that a total of 73 spectra distributed over the 5?3 x 4* region are 

sufficient to give a convolved spectrum, that is accurate to 

5 X 10^® cm"2 In integrated column density (see 5Ae). With this large number 
of spectra, the signal-to-nolse (S/N) ratio in is always very high 

even if only a little time is spent at each position. It thus makes sense to 

split the observations into two groups: target spectra , which are high S/N 

ratio data taken in the direction of Interest before and after map spectra 
which are taken by rapidly acannlng the area of the BL main beam closest to 
the direction of interest. The map spectra determine <S>, using equation Al. 
The target spectra are the onea which are actually corrected via equation A4. 

In section e) we discuss measurements of the rate of change of the stray 
radiation spectrum with time and position. The results indicate that target 
spectra can be observed for 30’’ to 40'*‘ before and after a map is made, and 
that it is generally desirable to center the map on a BL spectrum taken within 
1* of the target positions. 

d) Typical Stray Radiation Corrections for the 43 m 

The 43 m spectra from Figure 9 are ahown individually in Figure 10 along 
with the amount of each that is attributable to stray radiation (shaded 
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portions). There is so little HI In the main beam of the 43 m In this 
direction that the stray radiation Is a particularly large fraction of the 
signal: about 57Z of the observed upper spectrum and 48X of the lower 

spectrum can be attributed to atray radiation. Recall that both spectra are 
from the same direction, but were taken at different hour angles. 

Observations at a number of high latitude directions were corrected for 
stray radiation. The Integrated values of <S>, expressed as an equivalent HI 
column density, are plotted against the local sidereal time (LST) of 
observation In Figure 11. There Is a large scatter In the amount of stray 

radiation, but It reaches a minimum around 12^, when the galactic plane Is at 

low elevations. The amount of stray radiation Is not correlated with the 
total amount of hydrogen In the main beam. 

Figure 11 shows that several spectra apparently have ’’negative'* stray 
radiation, l.e., a signal that must be added to the 43 m spectra to bring them 

Into line with the BL spectra. The effect occurs only when the signal Is 

relatively large, and amounts to, at most, only 6Z of the total column 
density. It probably arises from a combination of baseline and calibration 
differences. The 43 m data were calibrated by adopting the Williams (1973) 
values for the brightness of the standard regions S6 and S8. The SBLH 
temperatures were calibrated against a cold load and appear to be a few 
percent higher than the Williams values. When the stray radiation In a 43 m 
map Is a negligible fraction of the total observed signal, either because the 
main beam emission Is quite bright or because of the LST and attitude of the 
telescope, the quantity <S> In equation A1 contains only the calibration and 
baseline differences between the 43 n and SBLH observations; It can easily be 
negative. Further Investigation of this phenomenon Is In progress, but It 
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doe* not appear to cause a aerious uncertainty In any of our data beyond the 

few percent level. 

e) T**t* of the Basic Aaeuraptione 

The tiro principal assumptions on which our correction procedure rests are 
(1) that the stray radiation spectrum does not change significantly over the 
time neceseary to make the 43 m observations, and (2) that the stray radiation 
in a single 43 m spectrum is well approximated by the <S> determined over 20 
square degrees from the convolved 43 m map (eq. A4). These assumptions, and 
other observations designed to test the accuracy of our results, are dlsct«saed 
below. 

1) Variation with Time 

For the first test about 75 hlgh-latltude directions were observed twice. 
In each case the second spectrum was taken about 40 min after the first. 
Differences between the spectra in a pair arise from noise, random baseline 
errors, and changes in the stray radiation component with time. The rma 
difference over the 75 pairs is only 4.0 x 10^^ cm~^ to which noise and 
baseline errors contribute typically < 2 x 10^^ cm'^. The shape of the 
spectra did not vary significantly over this period. 

As a second test, we produced three estimates of <Ti{{^q> from consecutive 
maps centered on almost the same position. All maps were made from a set of 
right ascension strips separated by 30' in declination. Data were averaged 
over 20* of time while the telescope was being driven la right ascension at 
either 90* /min (fast rate) or 30* /min (slow rate). A "slow" scan gives 
spectra e'^ery 10* in right ascension. A map et the fast rate was followed by 
one at the slow rate, offset 10* in declination, and then by another map at 
fast rate offset another 10* in declination. The entire procedure took about 
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two hours. Differences between the three resulting estimates of <T^ra(^ will 
be dominated by small angular scale structure in the HI, by time variations in 
<S>, and by baseline uncertainties. These observations were made in 15 
different directions. The maximum difference among the three values of 
<Tjiha(^ for a given direction averaged only 5.2 x 10^® cm“^, and was never 
larger than 1.2 x 10^^ cm"2. We conclude that <S> needs to be redetermined 
only evary two hours or so, and that maps made at the fast rate, which contain 
Just over 70 spectra, give an <T)^rao^ that is adequate for estimating the 
stray component in a 43 m spectrum. 
ii) Variation with Position 

The assumption that <S> derived over the -‘3* x 2* area of the BL main beam 
is representative of the stray radiation over any smaller region within that 
beam is reasonable, for far sidelobes are always very broad (see, e.g., KMR). 

We examined the variation in <S> that resulted from sequential 
obs'irvations of positions separated by 2*. Typical changes were 
< 4 X 10^® cm"2, although a change of 2.6 x 10^^ cm"2 was observed once. Over 
displacements of 4* to 10* the average difference in <S> wae 1.2 x 10^^ cm“2. 
These differences include the effects of noise and baseline uncertainties as 
well as the changes in <S> over times ^1 hour. 

Because these relatively large changes in position produce only modest 
changes in <S>, we conclude that angular structure in <S> on scales < 2* is 
modest, and will introduce only minor uncertainties into our results. 
ill) Internal Consistency 

The most useful estimate of the internal consistency of the stray 
radiation correction process comes from maps of The Hole, where pairs of 
43 m spectra, taken at sky positions separated by < 5* , but at times separated 


27 


by up to ••vcral d«yi, woro corrocttd uilni; different BL epectra. Ibe me 
deviation of the difference in the corrected for S3 auch pairs is 
7 X 10^^ ca'^t giving an average error, per spectnm, of 5 x 10^^ ca~^. This 
error includes noise and baseline uncertainties, and is at nost only lOZ of 
the total coluan density. At high latitude positions raore typical than The 
Role, the error would be < IZ of the total 

In addition, the stray radiation correction procedure was applied to data 
froa four observing sessions. TWo of th«; sessions, during March and August 
1983, used a receiver %rith dual linearly polarized scalar feeds. The other 
two sessions, i*> September 1983 and April 1984, used a receiver with dual 
linearly polarized hybrid-mode feeds. There was no significant polarization 
structure in the observed HI spectra, so data froa the two polarizations were 
always combined. Many directions were observed in three sessions, and three 
directions in all four. These observations give a check on the consistency of 
the correction technique over time and %d.th different antenna illumination 
patterns (and hence different main beaa efficiency and sidelobes). Table A1 
shows the average observed and corrected 43 a column densities for ten 
directions. The dispersion (lo) on each * measure of the total 

variation from all sources, including baselines and noise, over the sessions. 
The average dispersion in corrected spectra is 0.8 x 10^^ cm~^ or < 4Z of the 
average Ngx* average corrected a is close to the dispersion expected froa 

noise and baseline uncertainties. For the average, and eight of the ten 
indi idual directi(>ns, < Ogf,, indicating that the correction procedure 

has removed stray component whose magnitude fluctuates froa session to 
session. The stray radiation correction procedure thus gives spectra that are 
repeatable to a high degree of accuracy over several observing sessions and 
despite major equipment changes. 
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We conclude that the average error Introduced by- the assumption thaf^the — 
stray radiation spectrum Is constant across 2-3 degrees In the sky and over a 
period of up to two hours Is no greater than the 5 x 10^® cm”^ error In our 
data caused by noise and baseline errors. The maximum residual stray 
radiation error seems unlikely to be more than two or three times this value 
under any circumstance. 

Iv) Agreement with Other Determinations of Nht 

As discussed above, Kalberla, Mebold, and Reich (1980, hereafter KMR) have 
published HI profiles from the 100 m Effelsberg telescope that were corrected 
for stray radiation by a complete deconvolution of a model beam pattern using 
all-sky HI surveys. We observed two of the KMR directions; the corrected 43 m 
spectra at 21' resolution are compared with KMR spectra at 9* resolution In 
Figure 12. Toward KMR position 4 we find » 7.3 x 10^^ cm”2 and KMR find 
7.1 X 10l9 cm"2^ In very good agreement, although the KMR spectrum has a 
slightly higher peak temperature. Toward the north gr lactic pole (KMR 17) we 
find 9.2 X 10^^ cm"2 while KMR find only 6.7 x 10^^ cm“2. Most of the 
disagreement arises from the obvious baseline problems In the KMR spectrum; 
the profile shapes are In good agreement. Half or more of the observed signal 
has been removed from all four spectra In the stray radiation correction 
processes. 

Additional support for our method comes from comparison of as derived 
from La observations of distant hlgh-latltude stars with the corrected Njfj 
derived from 21 cm observations In the same direction. The six stars In the 
Lockroan, Hobbs, and Shull (1985) survey that are expected to lie entirely 
beyond the galactic HI layer all show ■ N2x to within the observational 
uncertainties In ('^ 15Z). This agreement provides Independent evidence 
that the overall 21 cm correction process gives accurate results. 
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f) Summary . ; .>iini...arv 

This summary Is based on our experience with the correction of x 10^ 
high latitude HI spectra for stray radiation. In general, corrected HI column 
densities Integrated over an entire spectrum should be accurate to 

X 10l9 cm~^ (lo) when all sources of error are Included. Higher accuracy, 
by a factor of ** 2, can be obtained In certain directions If special care Is 
taken In the data acquisition and especially In baseline removal (see Appendix 
B). Also, because most stray radiation comes from relatively low velocity gas 
that Is Doppler shifted by the combined earth-sun motion, spectral features at 
|v| > 100 km s~^ will usually not be contaminated and the accuracy of these Is 
limited only by Instrumental effects. We are continuing to study the stray 
radiation problem to refine the technique discussed here. 

The limiting factors on this technique probably reside In the SBLH 
spectra. These were taken with a velocity resolution of only 5.3 km s“^, so 
a narrow spectral line on the edge of a broad emission feature will appear 
shifted In velocity In their spectra. In extreme cases (of which we have a 
few examples). It is pointless to try to obtain corrected HI profiles at a 
velocity resolution higher than 5.3 km s“^. Another limitation Is that the 
SBLH survey covers only 6 > -AO*; more negative declinations cannot be 
corrected for stray radiation. Nonllnearltles In the SBLH filters (see SBLH) 
should not be a serious problem as these only affect bright lines In whose 
direction stray radiation Is usually a very small fraction of the total 
signal. Finally, our technique assumes that the SBLH spectra are totally 
correct. From examination of a number of these spectra we estimate that any 
one may have an uncertainty that Is equivalent to a column density error of 
2 X IQI® cm~^ ( lo) caused primarily by baseline errors. 
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Stray radlation' contrlbutlohB. to- th«'. SBLH»prof il-es-are estimated, to br i *T M " »-«•> t-f* 
typically leas than a few times 10^® cm”2, and at most 1 x 10^^ cm“2^ at 
Intermediate and high galactic latitudes. 

In summary, the method of "bootstrapping" high resolution spectra to the 
SBLH survey gives roughly an order of magnitude reduction in stray radiation 
over uncorrected high latitude HI spectra. This technique has the Important 
advantage over a formal deconvolution of the antenna response (as done by KMR) 
in that It can be easily applied to any telescope, with any receiver, since 
detailed knowledge of the full antenna pattern Is not required. 
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APPRENDIX B 

REDUCTION OF ERRORS IN BASELINE REMOVAL — - 

In the course of this project many thousands of HI spectra had to be 
corrected for Instrumental baseline* The baseline Is usually approximated by 
a low order polynomial which Is fit* In the least-squares sense, to the 
emission-free channels of each spectrum. Because the polynomial Is fit to 
noisy data, there are random errors In the derived coefficients which can be 
the dominant source of error In estimates of the power In a spectrum. 

This Appendix discusses a way to reduce the errors associated with baseline 
removal from sets of spectra. 

The area. A, under a spectral line. Is 

A - Av I (ti-Gi) B1 

L 

where t^ Is the observed temperature, Gj^ Is the value of the polynomial 
baseline at the 1^^ channel, Av Is the channel spacing (taken to be unity from 
now on), and the subscript L Indicates that the sum Is over the spectral 
channels which may ‘contain emission. If the polynomial G Is simply a constant 
go, then 

1 

H = So ” — I H . B2 

Nb B 

where Ng Is the number of channels that go Into the baseline fit and the 
subscript B Indicates that the sum Is over these channels. If all channels 
are Independent, and If each has the same variance then the variance of 

go Is given by 


(Ago)^ “ Ob^/Nb 


B3 
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In this case the area reduces to 

A - J ti - Nl go B4 

L 

where is the number of channels In the line, and the variance In the area 
follows from the law of propagation of errors: 


Aa 2 ■ 0^2 + og2 Nj^2/Ng 
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where a-i^ Is the variance of a channel containing line emission. Because 
OB^ • 0^2 and typically Ng ■ Nl, the uncertainty In the area has approximately 
equal contributions from noise In the signal and baseline regions of the 
spectrum. 

Observations often require fitting higher order baselines, such as the 
second order polynomials removed from our HI spectra. To estimate the 
associated errors we write the baseline as an orthogonal polynomial 

Gi - go + glB + g2 YiT2 + ••• B6 

where B ■ (1-B'), Yi - (1-Yi'), Y2 “ ~ Y2') B', Yl' and Y2* chosen 

to make gQ, g^ and g 2 orthogonal (e.g. Bevlngton 1969). The quantities B, Yl 
and Y 2 are functions only of the channel number and baseline regions. An 
extension of equation B5 gives the variance In the area: 


Aa2 




(I B)2 

L 


I b2 


+ 


(I YlY2)^ 

L 

+ ... 

I 

B 


B7 


All terms In equation B7 are positive, so typically the uncertainty In the 
area Is dominated by uncertainty In the baseline. 

This analysis suggests that the variance of the area can be reduced If 
more channels are Included In the baseline fit, l.e.. If both the bandwidth 
and total number of channels In a spectnim are Increased. This Is often not 
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feasible, and even when It is, wlde-band data generally need to be fit with « 
higher order polynomial baseline than narrow-band observations. 

It is more useful to recognize that in removing an instrumental baseline 
we are attempting to model a function that typically changes slowly with time, 
and is usually highly correlated between spectra taken sequentially at nearby 
positions. A baseline model obtained from a fit to the average of n, 
sequential spectra will have 03 ^ - oi^/rig and a more accurate Lk as long as 
the true baseline is well behaved. This information is not used if spectra are 
reduced individually. For second-order baselines an ng •• 3 is usually 
sufficient to make the terms of equation B7 containing 03 ^ smaller than the 
intrinsic uncertainty When ng > 1 the baseline that is removed from 

any spectrum is not optimum in the least-squares sense, but differences tend 
to be small and well within the range of uncertainty of the least-squares 
solution. 

The value of modeling the baseline from an average of ng spectra was 
tested using the HI data obtained for this paper. (Xir original observing 
procedure was such that we took 13 sequential spectra spaced over a 2* strip 
on the sky. The strip was repeated about 45''^ later, giving two spectra in 
each observed direction. The difference in area between the pairs of spectra 
gives a measure of the total uncertainty in the area from all sources. The 
average difference was attributed to a change in the stray radiation 
component. The dispersion arises from the baseline and noise. We reduced a 
hundred such pairs varying only the quantity ng between tests. An ng 1 is 
the usual procedure of fitting a baseline to an individual spectrum. For 
ng - 3 the baseline was fit to the average of three consecutive spectra and 
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then removed from the second spectrum of the triplet. • An ng > 5 was also . 
tested. As ng random baseline errors vanish and AA reaches Its lowest 
possible value. 

Table B1 shows the observed dispersion, AA, In K km s~^, compared with the 
theoretical dispersion as given by equation B7. The technique of averaging 
consecutive spectra to reduce the error In the baseline clearly works for 
these data, and It gives dispersions In reasonable agreement with theoretical 
expectations. Improvement In the quality of the baselines for ng > 1 Is 
especially noticeable on video displays of spectral line maps. 

This technique of spectrum averaging should be applicable to many 
different observations. The data presented In this paper were reduced with 
ng " 3, and while an ng « 5 did give some further reduction In AA, our 
conservative choice maximized the gain In accuracy while ralnlmlzlng the 
chances that a real baseline change occurred over the ng spectra. In all, some 
13000 spectra were reduced this way with satisfactory results. Some care must 
be taken In using a large ng with autocorrelation receiver spectra because 
these have ^tj^ ■■ 0 as a result of normalization. Observations of strong lines 
with a lot of spatial structure may require use of an ng no greater than 3 In 


these cases 
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TABLE 1 

OBSERVED DIRECTIONS- 
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45 

33 
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33 
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53 

59 

7.6 

182.0 

+58.4 

10 

27 

11 
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10 

53 
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12 

32 

44 

+28 

10 

45 
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. TABLE A1 


Consistency of Corrected Spectre 


*.b 
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Obs 

- Nhi ( 10-9 

• 

cra“2)— 

Cor 

• 


3 

38.1 1 

1.3 

35.1 

0.8 

66 -45 

3 

53.3 

1.0 

48.1 

0.4 

68 -63 

3 

29.5 

1.2 

24.4 

1.1 

88 -30 

3 

61.5 

3.6 

49.8 

0.1 

101 +65 

3 

15.5 

1.1 

12.1 

0.6 

182 +58 

4 

17.7 

1.2 

14.0 

0.6 

214 -32 

3 

63.1 

0.2 

55.6 

0.8 

237 -27 

4 

30.1 

1.9 

19.7 

0.9 

240 +37 

3 

39.0 

0.6 

38.0 

0.5 

316 +51 

4 

22.8 

1.6 

23.6 

2.3 
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TABLE B1 

Dispersion In an Area Estimate' Ar 


"s 

AA Thao. 
(K km s“l) 

/.A Obs. 

(K km 8-1) 

1 

3.5 

4.1 

3 

2.3 

2.4 

5 

2.0 

2.0 

m 
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FIGURE CAPTIONS 

Figure 1. The dlitributlon of the nuaber of spectre In e given column density 
bln for the feint end of the SBLH HI survey. Of the *^22000 spectra In that 
survey* only 0.9E have a total Nm < 1 x 10^0 cm"^, and < O.IZ have 
Nhi < 7 X 10^^ cn-2. These lowest column density directions (shaded portion 
of the figure) all come from a part of the sky near a ■■ 10^7* 6 ■ 58** that we 
call The Hole. 

Figure 2. The total HI column density toward The Hole. Positions are offsets 
from a - 10^ 45’’^* d * 57* 20*. -Column density contours are In units of 
1 X 1019 cm~^. Spectra were measured every 10* * In both coordinates, within 
the area outlined by the dashed border* The column density never drops 
significantly below 4.5 x 10l9 cm'^ over this map. 

Figure 3. Sample spectra toward The Hole showing the four line components. 
Offsets beside each spectrum refer to the scale In Figure 2. The spectrum at 
offsets (-65**0') Is the weakest galactic HI line detected to date. 


Figure 4. Column density maps for the Individual spectral components toward 
The Hole. Contours are In units of 1 x 10^9 cm~^* and a dashed contour 
Indicates the maximum (but uncertain) total extent of the emission. The 
positions are offsets from a ■ 10^ 45'** 6 ■ 57* 20* * as In Figure 2. The line 
components are defined In the text and Illustrated In Figure 3. Note the 
smooth variation In column density for the dominant low and Intermediate 
velocity cong>onentOf 

Figure 5. Line width (FVHM) and peak velocity maps of the cloud toward The 
Hole. Contours are In km a~^. 
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Pigur* 6* Th« reduced itructur* function, D* , vertuo angular separation, for 
the total N|fx, and for the Ngx in the low and Intemedlate velocity components 
toward The Hole. The figure shows the rms difference between column densities 
separated by 6, expressed as a fraction of the average column density In the 
component. It shows that nearby positions are highly correlated, and that the 
major structure Is In gradients. 

Figure 7. Column densities across eight low Ngj regions. Positions are 
offsets from the values given In Table 1. The dashed curve for the direction 
t ■ 182*, b " 58?4 shows the Ngx when the high velocity gas In this direction 
Is not included. There are no column densities < 7 x 10^^ cm*~2. 

Figure 8. Spectra in the low column density areas. These were taken at the 
positions listed in Table 1, and are at the center of each strip in Figure 7. 
Note the complexity of the typical spectrum. 

Figure 9. An example of the time variation of stray radiation. These 
spectra were taken at identical positions, but at two differert hour angles. 
The varying component Is only a fraction of the total stray radiation. 

Figure 10. The aame spectra as In Figure 9, along with the amount of each 
that arises from stray radiation (ahaded areas). About 57X of the upper 
spectrum and 48Z of the lower (pectrum are attributable to stray radiation. 

Figure 11. The amount of stray radiation, expressed as an equivalent column 
density of HI, that was removed from a number of high latitude observations, 

e 

plotted against the LST of observation. 
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Figure 12. Comparison of spectra at 9* angular resolution from the 100 m 
telescope corrected for stray radiation by an all-sky deconvolution (labeled 
KMR with the running number assigned in that paper) with spectra at 21' 
resolution obtained at the 43 m and corrected by the method described in this 
Appendix. Before correction for stray radiation the spectra contained about 
twice the signal shown here. Ihe agreement between the two methods is quite 
good, given the difference in angular resolution and the obvious baseline 
problems in the 100 m observations of KMR 17. 

















(8) so; 









cos(S) 






0+90 



Aa cos 8 



»00k Uf 0*90 









LSR 



OWQINAL PAGE IS 
OF POOR QUALITY 


8 









KMR 4 


Omo«NAL PAGE 
OF POOR QUALITY 


% 






